i 



MICROFICHE 
REFERENCE 
LIBRARY 

A project of Volunteers in Asia 

Building to Resist the Effect of Wind, Volume 2: 
Estimation of Ext reme Wind Speeds and Guide to 
the Det e rmination of Wind Forces 

by: Emil Simiu and Richard D. Marshall 

Published by: 

National Bureau of Standards 
U.S. Department of Commerce 
Washington, DC 20234 USA 

Paper copies are $ 1.30. Ask for stock number 
003-003-01718-3 when ordering. 

Available from: 

Superintendent of Documents 
US Government Documents 
Washington, DC 20402 USA 

Reproduction of this microfiche document in any 
form is subject to the same restrictions as those 
of the original document. 




NBS BUILDING SCIENCE SERIES 100 



Building To Resist The Effect Of Wind 

VOLUME 2. Estimation of Extreme Wind 

Speeds and Guide to the 
Determination of Wind Forces 



U.S. DEPARTMENT OF COMMERCE • NATIONAL BUREAU OF STANDARDS 



NBS BUILDING SCIENCE SERIES 100-2 



Building To Resist 
The Effect Of Wind 

In five volumes 

VOLUME 2: Estimation of Extreme Wind Speeds and 

Guide to the Determination of Wind 
Forces 

Emil Simiu 

Richard D. Marshall 

Center for Building Technology 

Institute for Applied Technology 

National Bureau of Standards 

Washington, D.C. 20234 



Sponsored by: 

The Office of Science and Technology 
Agency for International Development 
Department of State 
Washington, D.C. 20523 




U.S. DEPARTMENT OF COMMERCE, Juanita M. Kreps, Secretary 

NATIONAL BUREAU OF STANDARDS, Ernest Ambler, Acting Director 
Issued May 1977 



Library of Congress Catalog Card Number: 77-600013 
National Bureau of Standards Building Science Series 100-2 

Nat. Bur. Stand. (U.S.), BUg. Sri. Ser. 100-2, .29 pages (May 1977) 
CODEN: BSSNBV 



For sale by the Superintendent of Documents, U.S. Government Printing Office 
Washington, D.C. 20402 - Price $1.30 



Stock No. 003-003-01718-3 



ABSTRACT 



The Agency for International Development spon- 
sored with the National Bureau of Standards, a three 
and a half year research project to develop improved 
design criteria for low-rise buidings to better resist the 
effects of extreme winds. 

Project results are presented in five volumes. Volume 
1 gives a background of the research arti vities, ac- 
complishments, results, and recommendations. In 
Volume 3, a guide for improved use of masonry 
fasteners and timber connectors are discussed. 
Volume 4 furnishes a methodology to estimate and 
forecast housing needs at a regional level. Socio- 
economic and architectural considerations for the 
Philippines, Jamaica, and Bangladesh are presented in 
Volume 5. 

Volume 2 consists of two reports. The first reviews the 
theoretical and practical considerations that are perti- 
nent to the estimation of probabilistically defined 
wind speeds. Results of the statistical analysis of ex- 
treme wind data in the Philippines are presented and 
interpreted. Recommendations based on these results 
are made with regard to the possible redefinition of 
wind zones, and tentative conclusions are drawn 
regarding the adequacy of design wind speeds cur- 
rently used in the Philippines. Report two describes 
some of the more common flow mechanisms which 
create wind pressures on low-rise buildings and the 
effects of building geometry on these pressures. It is 
assumed that the basic wind speeds are known and a 
procedure is outlined for calculating design wind 
speeds which incorporates the expected life of the 
structure, the mean recurrence interval, and the wind 
speed averaging time. Pressure coefficients are tabu- 
lated for various height-to-width ratios and roof 
slopes. The steps required to calculate pressures and 
total drag and uplift forces are summarized and an il- 
lustrative example is presented. 



Key words: Building codes; buildings; codes and standards; 
housing; hurricanes; pressure coefficients; probability dis- 
tribution functions; risk; statistical analysis; storms; struc- 
tural engineering; tropical storms; wind loads; wind speeds. 
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1. ESTIMATION OF EXTREME 
WIND SPEEDS- 
APPLICATION TO THE 
PHILIPPINES 

by 

E. Simiu 
1.1 INTRODUCTION 

In modern building aides and standards | 1 . 2| basic 
design wind speeds are specified in explicitly pro- 
babilistic terms. At any given station a random varia- 
ble can be defined, which consists of the largest yearly 



wind speed. It the station is one for which wind 
records over a number of consecutive wars are 
available, then the cumulative disti ibution (unction 
(CI)[)o; this random \ariable maw.it least in theory 



be estimated to characterize the probabilistic behavior 
of the largest yearly wind speeds. The basic design 
wind speed is then defined as the speed correspond- 
ing to a specified value Fq of the CDF or, equivalently 
(in view of the relation N= 1/ (1-Fq) in which N = 
mean recurrence interval), as the speed correspond- 
ing to a specified mean recurrence interval. For exam- 
ple, the American National Standard A58.1 [1] 
specifies that a basic design wind speed corresponding 
to a 50-year mean recurrence interval (i.e , io a value 
Fq of the CDF equal to 0.98, or to a probability of ex- 
ceedance of the basic wind speed in any one year 
equal to 0.02) be used in designing all permanent 
structures, except those structures with an unusually 
high degree of hazard to life and property in case of 
failure, for which a 100-year mean recurrence inter- 
val (Fq = 0.99) must be used, and structures having no 
human occupants or where there is negligible risk to 
human life, for which a 25-year mean recurrence 
(Fq = 0.96) may be used. A wind speed corresponding 
to a N-year recurrence interval is commonly referred 
to as the N-year wind. 

The mean recurrence intervals specified by building 
codes, rather than being based on a formal risk 
analysis — which is in practice not feasible in the pre- 
sent state of the art — are selected in such a manner as 
to yield basic wind speeds which, by professional con- 
sensus, are judged to be adequate from a structural 
safety viewpoint. Nevertheless, it is generally 
araimed that adequate probabilistic definitions of 
design wind speeds offer, at least in theory, the ad- 
vantage of insuring a certain degree of consistency 
with regard to the effect of the wind loads upon struc- 
tural safety. This is true in the sense that, all relevant 
factors being equal, if appropriate mean recurrence 
intervals are used in design, the probabilities of failure 
of buildings in different wind climates will, on the 
average, be the same. 

In the practical application of the probabilistic ap- 
proach to the definition of design wind speeds, cer- 
tain important questions arise. One such question per- 
tains to the type of probability distribution best suited 
for modeling the probabilistic behavior of the extreme 
winds. The provisions of the National Building Code 
of Canada [2j are based upon the assumption that this 
behavior is best modeled by a Type I (Gumbel) dis- 
tribution. The American National Standard A58.1 [1], 
on the other hand, assumes that the appropriate 
models are Type II (Frechet) distributions with loca- 
tion parameters equal to zero and with tail length 
parameters dependent only upon type of storm. 
Finally, Thorn 129] has proposed a model consisting of 
a mixed probability distribution, the parameters of 
which are functions of (a) the frequency of occurrence 
of tropical cyclones in the 5° longitude-latitude square 
under consideration and (b) the maximum average 



monthly wind speed recorded at the station investi- 
gated. The question of selecting the most appropriate 
distribution is one that deserves close attention: in- 
deed, as indicated in References 23 and 22, the mag- 
nitude of the basic design wind speed may depend 
strongly upon the probabilistic model used. 

Assuming that the type of probability distribution best 
suited for modeling the behavior of the extreme 
winds is known, a second important question arises, 
viz., that of the errors associated with the probabilistic 
approach to the definition of design wind speeds. 
Such errors depend primarily upon the quality of the 
data and upon the length of the record (i.e., the sam- 
ple size) available for analysis. 

These questions will be dealt with in this work, which 
will also present results of statistical analyses of wind 
speed data recorded in the Philippines. In the light of 
the material presented herein, possible approaches 
will be examined to the definition of extreme wind 
speeds for purposes of structural design in the Philip- 
pines. 

1.2 WIND SPEED DATA 

For the statistical analysis of extreme wind speeds to 
be meaningful, the data used in the analysis must be 
reliable and must constitute an homogeneous set. The 
data may be considered to be reliable if: 

• The performance of the instrumentation used 
for obtaining the data (i.e., the sensor and the record- 
ing system) can be determined to have been adequate. 

• The sensor was exposed in such a way that it 
was not influenced by local flow variations due to the 
proximity of an obstruction (e.g., building top, ridge 
or instrument support). 

A set of wind speed data is referred to herein as 
homogeneous if all the data belonging to the set may 
be considered to have been obtained under identical 
or equivalent conditions. These conditions are deter- 
mined by the following factors, which will be briefly 
discussed below: 

• type of instrumentation used 

• averaging time (i.e, whether highest gust, fastest 
mile, one-minute average, five-minute average, etc. 
was recorded). 

• height above ground 

• roughness of surrounding terrain (exposure) 

• in the case of tropical cyclone winds, distance 
inland from the coastline. 

1.2.1 Type of instrumentation 

If, during the period of record, more than one type of 
instrument has been employed for obtaining the data, 
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the various instrument characteristics (anemometer 
and recorder) must be carefully taken into account 
and the data must be adjusted accordingly. 

1.&2 Averaging Time 

If various averaging times have been used during the 
period of record, the data must be adjusted to a com- 
mon averaging time. This can be done using graphs 
such as those presented in Reference 19 and included 
in figure 1 in which Z is a parameter defining the 
terrain roughness (see, for example, Ref. 10). 

1.2.3 Height Above Ground 

If, during the period of record, the elevation of the 
anemometer had been changed, the data must be ad- 
justed to a common elevation as follows: Let the 
roughness length and the zero plane displacement be 
denoted by Z 8 and Z d , respectively (Z„, Z d , are 
parameters which define the roughness of terrain, see 
Ref. 10). The relation between the mean wind speeds 
IMZJ and U(ZJ over horizontal terrain of uniform 
roughness at elevation Z t and Z % above ground, 
respectively, can be written as 



U(Z.) 



En z j~ z d 




(1) 



Suggested values of the roughness length Z„ are given 
in table 1 (see refs. 10,21,7). For example, at Sale, 
Australia, for terrain described as open grassland with 
few trees, at Cardington, England, for open farmland 
broken by a few trees and hedge rows, and at 
Heathrow Airport in London, Z = 0.08 m [ 10, 21 ]. At 
Cranfield, England where the ground upwind of the 
anemometer is open for a distance of half a mile 
across the corner of an airfield, and where neighbor- 
ing land is broken by small hedged fields, Z„ = 0.095m 
[91. The values of Z„ for built-up terrain should be 
regarded as tentative. It is noted that Equation 1 is ap- 
plicable to mean winds and should not be used to 
represent the profiles of peak gusts. 

Table 1. Suggested Values of Z for Various Types 
of Exposure 



Type of Exposure 


Z (meters) 


Coastal 


0.005-0.01 


Open 


0.03- 0.10 


Outskirts of towns, suburbs 


0.20- 0.30 


Centers of towns 


0.35- 0.45 


Centers of large cities 


0.60- 0.80 



The zero plane displacement Z d may in all cases be 
assumed to be zero, except that in cities (or in wooded 
terrain) Z d = 0.75 h, where h = average height of 
buildings in the surrounding area (or of trees) [10, 16]. 
Thus, for example, if in open terrain with Z = 0.05 m, 
LH.23) = 30 m/s, then adjustment of this value to the 
height Z, = 10 m, using Equation 1, gives 



UUO) = U(23) 



0. 05 
0. 05 



= 30 



5.30 
6.13 



25.9 m/s. 



It is noted that, in most cases, the roughness 
parameters Z , Z d must be estimated subjectively, 
rather than being determined from measurements. 
Good judgment and experience are required to keep 
the errors inherent in such estimates within -easona- 
ble bounds. In conducting statistical studies of the ex- 
treme winds, it is advisable that for any particular set 
of data, an analysis be made of the sensitivity of the 
results to possible errors in the estimation of Z and 

In the case of winds associated with large-scale ex- 
tratropical storms, the mean wind LKZ) at height Z in 
terrain of roughness Z , Z rf is related as follows to the 
mean wind l/,(Z,) at height Z in terrain of roughness 
Z 01 ,Z dl [2ll: 



U(Z) = j3 U.(Z.) 



(2) 



fin 



Z.-Z 



The quantity, (J, may be obtained from figure 2, 
which was developed in Reference 21 on the basis of 
theoretical and experimental work reported by 
Csanady [4] and others [26]. (Note that Z 0I <Z .) 

Equation 2 may be applied if the roughness of the ter- 
rain is homogeneous over a horizontal distance from 
the anemometers of about 100 times the anemometer 
elevation [18, 24]. 

Let, for example, U(Z,) = 29 m/s, Z- = 10 m, Z^ = 0.05 
m, Z dl = 0. The corresponding speed U(Z) at Z = 40 m, 
say, in open terrain of roughness Z> = 0.25 m, Z d = is 



fin — 

U(10) = 1.12 X 29— Pjp =31.1 m/s. 

fin — 

0.05 



where 1.12 is the value of for Z 01 = 0.05 m, Z„ = 0.25 
m, obtained from figure 2. 
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It is pc "nted out that, just as in the case of Equation 1, 
errors are inherent in Equation 2 that are associated 
with the subjective estin u> , .ion of the roughness 
parameters. Also, recent research suggests that in the 
case of tropical cyclone winds Equation 2 underesti- 
mates wind speeds over built-up terrain, calculated as 
functions of speeds over open terrain, by amounts of 
the order of 15% or more[17]. 

1.2.4 Distance Inland from the Coastline 

The intensity of hurricane or typhoon winds is a 
decreasing function of the distance inland from the 
coastline. Hurricane wind speeds may be adjusted to a 
common distance from the coastline by applying 
suitable reduction factors. Such reduction factors have 
been proposed by Malkin, according tc whom the 
ratios of peak gusts at 48, 96 and 144 km from the 
coastline to peak gusts at the coastline are 0.88, 0.82 
and 0.78, respectively [8, 14]. 

1.3 PROBABILISTIC MODELS OF EXTREME 
WIND SPEEDS 

The nature of the variate suggests that an appropriate 
model of extreme wind behavior is provided by prob- 
ability distributions of the largest values, the general 
expression for which is 111]: 

F(f) = exp | -\(v-li)lvl~ y f M < » < °° 

-oo<Ji<oo (3) 
0<O"<oo 

y>0 

where v = wind speed, (i — location parameter, <r — 
scale parameter, y = tail length parameter. Equation 3 
may be regarded as representing a family of distribu- 
tions, each characterized by a value of the tail length 
parameter y. As y becomes larger, the tail of the pro- 
bability curve becomes shorter, i.e., the probability of 
occurrence of large values of the variate becomes 
smaller. In particular, as y — ► oo. Equation 3 may be 
shown to become 

F(i>) = exp | -exp[-{v-fi)l<r\ | -oo <u< oo 

— oo <n < oo (4) 
0<<r<oo 

The distributions given by Equations 3 and 4 are 
known as the type II and the type I distributions of 
the largest values, respectively. 

Two basic procedures for estimating probabilities of 
occurrence of extreme winds are currently in use. The 
first procedure consists in estimating the parameters 
of a probability distribution of the largest values from 
the series of annual highest wind speeds at the station 
considered. This procedure has been applied by 



various authors as follows: 

(a) In Reference 23, estimates are made of all three 
parameters, fi, a- and y in Equation 3, no specific 
value being assigned a priori to any of these 
parameters. 

(b) In References 27 and 28, the location parameter is 
assumed to have the value n = 0. Estimates are then 
made of the remaining parameters, o- and fi. The ar- 
bitrary assumption that n = entails a sacrifice in 
goodness of fit; the justification for using this assump- 
tion is that it makes possible the application of 
Lieblein';. well-known estimation procedure [13] to 
obtain values of <r and y [27]. However, in view of the 
recent development of alternative estimation pro- 
cedures applicable to type II distributions with n*0 
[23], the assumption that fi = becomes unneces- 
sary. 

(c) Court in the United States [3], Davenport in 
Canada [5] and Kintanar in the Philippines [12] have 
assumed the universality of the type I distribution, 
i.e., that the tail length parameter is always y = oo. 
Estimates are then made of the parameters fi and o\ 

The second procedure assumes the universal validity 
of the mixed distribution 

Hv) = p £ exp [-£) " '] + p T expH£) (5) 

proposed by Thorn in Reference 29. The first and the 
second term in the sum of Equation 5 represent the 
probabilities that the winds associated with extra- 
tropical storms and with tropical cyclones, respec- 
tively, will not exceed the value, v, in any one year. 
The scale parameter, <r, is an explicit function of the 
maximum of the average monthly wind speeds 
recorded at the station considered. The second 
parameter of the mixed distribution, pj,\s an explicit 
function of the frequency of occurence of tropical 
cyclones in the 5° longitude-latitude square under con- 
sideration, and p E = 1 - pj. Thus, in this second pro- 
cedure, the series of annual highest winds is not used 
for estimating distribution parameters. 

An assessment of the models described in this section 
will now be presented. 

1.4 ASSESSMENT OF PROCEDURES BASED 
ON THE ANNUAL HIGHEST SPEEDS 

To assess the validity of current probabilistic models, 
statistical analyses of annual highest speeds were car- 
ried out using a computer program described in 
Reference 23. The results of the analyses, which are 
reported in detail in Reference 23, lend credence to 
the belief that a sufficiently long record of annual 
largest speeds will provide an acceptable basis for 
probabilistic extimatesof the N-year winds — even for 
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large values of N, such as are of interest in structural 
safety calculations — if the following conditions are 
satisfied. First, the value of opt {y) for that record is 
large, say y > 40 (opt(y) = value of y [see eq. 3] for 
which the best distribution fit of the largest values is 
obtained). Second, meteorological information ob- 
tained at the station in question, as well as at nearby 
stations at which the wind climate is similar, indi- 
cates that winds considerably in excess of those 
reflected in the record cannot be expected to occur ex- 
cept at intervals many times larger than the record 
length. Wind climates which satisfy these two condi- 
tions will be referred to as well-behaved. 

Assuming that the wind speed data are reliable, lower 
bounds for the sampling error in the estimation of the 
N-year winds in a well-behaved climate may be 
calculated on the basis of a mathematical result, the 
Cramer-Rao relation, which states that for the type 1 
distribution (see ref. 11, p. 282) 

1.10867 , 
var(/i) > a- 2 

n (6) 

... 0.60793 , 
Var (<r) > or 2 

n (7) 

where var (£), var («r) are the variances of p. and o-, 
where p. and & are the estimated values of n and <r, 
respectively, obtained by using any appropriate 
estimator consistent with basic statistical theory re- 
quirements; er is the actual value of the scale 
parameter and n is the sample size. Using Equations 6 
and 7, lower bounds for the standard deviation of the 
sampling error in the estimation of the N-year wind, 
SEflp(N)], can be approximated as follows. Equation 4, 
in which the parameters (t, <r are replaced by their 
estimates n, a- , is inverted to yield 

v(N) = /i-ai-Jj)3- (8) 

where 

G(l-Jj) = -ln[-ln(l-£j)]= (9) 

Then 

SD[v(N)] > [varifi) + [G(l-^)] 2 var (*)]' /2 (10) 

Equation 10 is based on the assumption that the error 
involved in neglecting the correlation between fi and 
<r is small. The validity of this assumption was 
verified by using Monte Carlo simulation techniques. 

Since the actual value of o- is not known, in practical 
calculations the estimated value cr is used in Equations 
6 and 7, For example, the distribution parameters cor- 
responding to the wind speed data at Davao (n = 24, 



see table 2), estimated by using the technique 
described in Reference 23, are p. = 38.89 km/hr, & = 
9.40 km/hr. It follows from Equation 8 that v (50) = 75 
km/hr and from Equations 6, 7, and 10 that SD[z> (50)] 
_> 5.18 km/hr. Subsidiary calculations not reported 
here have shown that Equation 10 provides a good 
indication of the order of magnitude of the sampling 
errors. 

1.4.1 Wind Climates Characterized by Small 
Values of opt (y). 

Occasionally, a record obtained in well-behaved 
wind climates may exhibit small values of opt (-y); this 
will occur if that record contains a wind speed that 
corresponds to a large mean recurrence interval. 
There are regions, however, in which, as a rule, the 
statistical analysis of extreme wind records taken at 
any one station yields small values of opt (y). This is 
the case if, in the region considered, winds occur that 
are meteorologically distinct from, and considerably 
stronger than the usual annual extremes. Thus, in the 
regions where tropical cyclones occur, opt (y) will in 
general be small, unless most annual extremes are 
associated with tropical cyclone winds. An example 
of a record for which y (opt) is small is given in figure 
3a, which represents the probability plot with y = opt 
(y) = 2 for the annual extreme fastest mile-speeds 
recorded in 1949-73 at the Corpus Christi, Texas, air- 
port. For purposes of comparison, the same data have 
been fitted to a type I distribution (opt ( y) = oo, or Eq. 
4); the fit in this case is seen to be exceedingly poor, 
i.e., the plot deviates strongly from a straight line (fig. 
3b). As shown in Reference 23, a measure of the good- 
ness of fit is given by the extent to which the pro- 
bability plot correlation coefficient is close to unity; 
this coefficient is printed out in figures 3a and 3b. 

To small values of the tail length parameter there fre- 
quently correspond implausibly high values of the 
estimated speeds for large recurrence intervals. In the 
case of the 1912-48 record at Corpus Christi, for exam- 
ple, opt (y) = 2 and the estimated 5-minute average is 
327 mph (155 m/s) for a 1000-year wind, which is 
highly unlikely on meteorological grounds. For 20- 
year records, the situation may be even worse: thus, 
for the 1917-36 Corpus Christi record, which contains 
an exceptionally high wind speed due to the 1919 
hurricane 13, 25], opt (y) = 1 and the calculated 1000- 
year wind is 1952 mph (873 m/s) [23], a ridiculous 
result. Also, the situation is not likely to improve sig- 
nificantly if the record length increases. From a 74- 
year record, a plot quite similar to figure 3 would 
presumably be obtained, with twice as many points 
similarly dispersed, to which there would correspond 
a similar least squares line on probability paper. 

It may be stated, consequently, that while in the case 
of well-behaved climates it appears reasonable to in- 
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fer from a good fit of the probability curve to the data 
that the tail of the curve adequately describes the ex- 
treme winds, such an inference is not always justified 
if opt (y) is small. 

It may be argued that one could avoid obtaining 
unreasonable extreme values by postulating that the 
annual largest winds are described by a probability 
distribution of the type I, i.e., by assigning the value y 
= oo to the tail length parameter. This has been done 
by Court l3l and Kintanar [12]. As can be seen in 
figure 4, the corresponding fit may be quite poor. 
However, the estimated extremes at the distribution 
tails will be reduced. The drawback of this approach 
is that unreasonably low estimated extremes may be 
obtained. For example, at Key West, Florida, if all 
three parameters of Equation 3 are estimated as in 
Reference 23, to the 1912-48 record there corresponds 
v (100) = 99 mph (44.2m/s) and v (1000) = 188 mph (84 
m/s)- see Reference 23. If it is postulated that y = oo, 
then v(50) = 70 mph (31.7 m/s), 1/(100) = 77 mph (34.4 
m/s) and u(1000) = 97 mph (38.8 m) [23l, an unlikely 
result in view of the high frequency of occurrence of 
hurricanes (about 1 in 7 years) at Key West. 

It may also be argued that since the estimated ex- 
tremes resulting from small values of y (say y < 4) 
may be too large, and those corresponding to y = oo 
may be too small, a probability distribution that might 
yield reasonable results is one in which y has an in- 
termediate value, say 4< y < 9. Such an approach 
has been proposed by Thorn and will now be ex- 
amined. 

1.5 ASSESSMENT OF PROCEDURE BASED 
ON THE HIGHEST AVERAGE 
MONTHLY SPEED 

The procedure for estimating extreme winds in hur- 
ricane-prone regions on the basis of annual highest 
winds at a station was seen to have the following 
shortcomings. First, because hurricane winds are 
relatively rare events, the available data may not con- 
tain wind speeds associated with major hurricane oc- 
currences and are therefore not representative of the 
wind climate at the station considered (see the case of 
Calapan in Section 1.7 of this report). Second, in 
regions subjected to winds that are meteorologically 
distinct from, and considerably stronger than the 
usual annual extremes, implausible estimates may be 
obtained. 

The model proposed by Thorn lEq. 5] in Reference 29 
represents an attempt to eliminate these shortcom- 
ings. It can be easily shown by applying the inter- 
mediate value theorem, that if this model is assumed, 
the estimated extreme winds may be obtained by in- 
verting an expression of the form: 

F(»)=expl(-v/o-)-7<t')] 



in which 4.5 < y (v) < 9. If the mean rate of arrival of 
tropical cyclones in the region considered is high, 
then y ( p) will be closer to 4.5. Otherwise, yiv) will be 
closer to 9; in regions where hurricanes cannot he ex- 
pected to occur, y (v) = 9. In order that estimates not 
be based upon possibly unrepresentative annual ex- 
treme data. Thorn's model does not make use of an- 
nual extreme speeds. Rather, the parameter cr is esti- 
mated from the maximum of the average monthly 
wind speeds on record at the location considered, 
presumably a quantity for which the variability is 
small. 

While the quasi-universal elimatological distribution 
proposed by Thorn is tentative, it will yield results 
which, for a first approximation, may in certain cases 
be regarded as acceptable. This model has recently 
been used by Evans l6j as a basis for obtaining design 
wind speeds for Jamaica. Estimates of extreme speeds 
obtained by Evans are substantially higher than the 
results obtained by Shellard [20] in his 1971 analysis 
of Commonwealth Caribbean wind data. 

It was shown in the preceding section that the ap- 
proach which utilizes the series of annual largest 
speeds may fail in regions in which hurricanes occur. 
For such regions, therefore, it may be that alternative 
approaches need to be developed. Among such ap- 
proaches is one in which estimates of extreme winds 
are based upon the following information: 

• average number of hurricanes affecting the 
coastal sector considered (per year) 

• probability distribution of hurricane intensities 

• radial dimensions of hurricanes 

• dependence of wind speeds upon central 
pressure and distance from hurricane center. 

This approach appears to provide useful estimates of 
extreme winds corresponding to large recurrence in- 
tervals — which are of interest in ultimate strength 
calculations — and is currently under study at the Na- 
tional Bureau of Standards. 

1.6 STATISTICAL ANALYSIS OF EXTREME 
WIND DATA IN THE PHILIPPINES 

Through the courtesy of the Philippine Atmospheric, 
Geophysical and Astronomical Services Administra- 
tion (PAG ASA), 16 sets of data were obtained consist- 
ing of maximum yearly wind speeds recorded during 
at least 14 consecutive years. The data for each of the 
16 stations are listed in table 2. Table 3 includes sub- 
jective station descriptions provided by PAG ASA 
personnel and the results of the analysis. In Table 3 
are listed VrfP^V^ W-year wind based on the dis- 
tribution for which the best fit of the largest values is 
obtained and V N °° = N-year wind based on the type I 
distribution, N = mean recurrence interval in years. 
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TABLE 2. MAXIMUM ANNUAL WINDS (ONE MINUTE AVERAGES) 



NO. 


Station 


Period of Record 


Maximum Annual Winds for Each Year of Record (km/hour) 


1 


Davao 


1950-73 


39,52,40,39,40,37,35,35,32,40,40.40,80,48,48,48,56,46,52,50,46,52,46,46 


2 


Cagayan de 
Oro 


1950-73 


47,24,19,13,19,19,12,12,12,19,16,14,21,6,24,17,19,37,37,46,37,48,41,41 


3 


Zamboanga 


1950-73 


in / ( mr\ 1t~\ AO f 1 At Af\ AO AC AO AO AO Cf\ C£. £-0 CI "7/1 C/l il 7^ £1 70 

48,64,40,3" 48,6 1 ,43,40,48,4 !^4oy2,4o,4e,DU,3o,OD,o7,/U,!ro,o 1 ,/ 4,o 1 ,/ o 


4 


Pasav Citv 


1950-73 


89,103,89,92,72,72,64,72,W,7z,ol,oo,oV,/4,l JU,b3,oU,l 1 1,oj,/4,^UU,oU,1 1 1,do 


5 


Manila 


1949-70 


■» nr r»"7 on i ni fV7 i nn i nc o i *7~\ o7 oo 1 1 1 1 fiC 1 fin 1 CO. *TA QQ 1 m 111 OH "if\f\ 

72,105,S 1 7,89,I0i,y7, 100,105,01,72,77,0*7, I21,103,100,lbo,/4,O7,lU/,l 1 l,yo,2UU 


6 


Manila 
Central 


1902-40 


46,56,65,80,73,55,77,70,41 ,68,50,69,64,68,42, 41 ,67,54,70,83,58,53,60,5 1 ,45, 

t -I cn TA M CC CO 11 1 t\f\ t(l CT 1 m C£ 

63,50,70,52, 55,58,37, 100,6U,45,52, 103,53, 36 


7 


Mirador 


1914-40 


79,79.7079,121,102,94,72,105,107,122,58,89,107,93,63,77,92,63,53,93,73,75, 
49 39 68 83 


g 


DclbU IU 


1950-73 


"97, 165,56,97,97,81, 64,64,6 1,48,97, 48,53,9(4,1 1 1,107,80,144,1(17,102,137,85,89,133 


9 


Calapan 


1959-73 


145,185,97,72,40,68,40,96,109,41,33,1 1 1,102,1 1 1,83 


10 


Surigao 


1952-73 


105,43,1 1 3,40,40,40,43,48,64,64,89,39,74,52,78, 104, 167,1 1 1,96,107,96,74 


11 


Laoag 


1949-73 


34,72, 1 1 8,71 , 1 08, 1 00,64, 8 1 ,8 1 ,64,8 1 ,79,64, 1 05,90, 1 44, 1 44,78, 1 20, 1 37, 1 (10,67 ,89,n7,70 


12 


lloilo 


1949-73 


97,106,64,64,97.64,64.64,72,74.64.7878,89,74, : 61 ,52,6 1 , 89,56,74 


13 


Cebu 


1950-73 


97,56,1 2 1 ,48,8 1 ,64,48,48,48,43,48,43,64,64,56,48,56,70,93,65,56, 10074,59 


14 


Legaspi _ 


1955-73 


40,97,97,48, 1 29,1 85,97,82,89, 1 04,74,89, 1 48,74,70, 1 74, 1 48.204,8 1 


15 


Tacloban 


1949-73 


137,58,106,1 13,56,68,60,47,48,69,87,42,105,93,74,1 II, 
70, 1 94, 1 33, 1 67,63,8 1 , 1 55. 1 04.67 


16 


Infanta 


1960-73 


60,43,45.50.128,39.133,126,46,46,189,85,104,83 



TABLE 3. STATION DESCRIPTIONS AND ESTIMATED EXTREME WIND SPEEDS 



No. 


SUIlan 


Wins 
/one — 

See Ret. 15 


Period of 
Record 


No. of 
Years 


Anemometer 
Elevation (meters) 


Description 
of Terrain 


opt (y) 


v" 


^^'•'(km/hrl 


V™' (km/hr) 


N=50 


N=100 


N = IUOU 


N=50 


N=10u 


N=1000 


I 


Davao 


III 


195073 


24 


Hi 


Town 










75 


S3 


las 


2 


Cagayan 
de Oro 


III 


1950-73 
















61 


69 


94 


3 


Zamboanga 


III 


1950-73 


24 


II) 


Open 


oo 








88 


95 


117 


4 


i'asav (-itv 


II 


1950-73 


24 


ll) 


Airport 




:30 


595 


A50 


■ I'M 


ISO 


534 


5 


Manila 


II 


l«M<a-7tl 


22 


Tup of 5-story 
building 


Port Area 


7 


2DH 


238 


3*5 


192 


212 


276 


6 


Manila 
Central 


ll 


1902-40 


39 


6 a 


Open (Park) 


35 


103 


112 


144 


101 


111) 


138 


I - 


Miradur 
Baguto 


11 


1114-411 


27 


16* 

(1500 m above 


Mountain top 
sea level) 










139 


151 


192 


8 


Baguio 


II 


195(1-73 


24 


1(1 


Muutain tup 
(150(1 m above 
sea level) 










164 


180 


230 


9 


Calapan 


II 


1959-73 


15 


iO 


Tup of hill; 
tuwn on une 
side. 










209 


234 


316 


ib - 


Surigao 


ll 


1952-73 


22 


7 


Large City 
Top of 40m 
hill 


40 


170 


191 


262 


167 


187 


250 


11 


laoag 


11 


1949-73 


25 


12 


Open 










174 


192 


252 


12 


lloilo 


II 


1950-73 


24 


25 (5m above 
20m bids.) 


Town 


2 


180 


235 


640 


■T38 


rsi - 


192 


13 


Cebu 


II 


1950-73 


24 


10 


Airport 


oo 






127 


140 


185 




|4 


ffiff 


1 


1955-73 


19 


12 


Airport*- 


90 


235 


264 


350 


-S5" 


264 


340 


T5— 


lacloban 


I 


1949-73 


■55 


12 (3m above 
9m bldg.) 


lopul 15m 
hill Military 
Camp" 1 


14 


213 


242 


352 


-204 


-BIT 


306 


I! - 


Intanta 


I 1 


19*0-73 


14 


10 


Town; Residen- 
tial area 


6 


T45 - 


290 


498 


214 


242 


533 



a 3 cup anemometer; mean speed averaged over one minute. 
''Mean speed averaged over one minute. 
c Trees at East side of anemometer. 



^North and East: sea exposure. 
e Omitted if opt (y) = oo. 
' One minute averages. 



1.7 INTERPRETATION OF RESULTS 

The results will be grouped into three classes, accord- 
ing to the wind zone (as defined in Ref. 15) in which 
the stations considered are located (table 3). 

1.7.1 Zone III 



It is noted that for all three Zone HI stations listed in 
Table 3, opt (y) = oo. It is convenient to adjust the 
speeds at Davao and Cagayan de Oro to open terrain 
exposure. On the basis of the terrain descriptions of 
Table 3, if it is assumed Z = 0.30 m, Z d = 0, Z 01 = 0.08 
m, Z^, = 0, it follows from Equation 2 that 
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M1°L . 8 
U(10) 

where U(10), U^IO) are mean speeds above ground in 
town and open exposure, respectively. Thus, in 
Davao and Cagayan de Oro the calculated 50-year 
mean speeds at 10 m above ground in open terrain are 
58 mph and 47 mph (94 km/hr and 76 km/hr) respec- 
tively, versus 55 mph, (88 km/hr) in Zamboanga. If the 
corresponding highest gusts are obtained by multiply- 
ing the one-minute means by a factor of, say, 1.20 (see 
fig. 1>, the estimated highest 50-yr gusts at Davao, 
Cagayan de Oro and Zamboanga are at most 94 x 1 .20 
= 113 km/hr, (70mph), i.e., considerably lower than 
the value specified for design purposes by the Na- 
tions' Structural Code of the Philippines [15] for Zone 
III, viz., 95 mph (153 km/hr). This suggests that the re- 
quirements of Reference 15 regarding wind loading in 
the Zone III portion of Mindanao are conservative 
and might be somewhat reduced. (It can be easily 
shown, on the basis of Eq. 2 and figure 1, that this 
statement holds even if it is assumed that the errors in 
the estimation of the parameter values Z = 0.30 m 
and Z = 0.08 m are of the order of as much as 50%.) 
To validate such a conclusion it would however be 
necessary to determine, from long-term records of 
tropica] cyclone occurences, that the 1950-73 data at 
the three stations analyzed are indeed representative 
for southern Mindanao. 

1.7.2 Zone II. 

Several difficulties arise in interpreting the results for 
the Zone II stations in table 3. It is noted, first, that the 
results obtained at stations in and near Manila (sta- 
tions 4, 5, 6 in table 3) are widely divergent. The dis- 
crepancies between the results for Pasay City and 
Manila may be due to the different elevations of the 
respective anemometers. It may also be conjectured 
that the discrepancies between these results and those 
obtained from the 1902-1940 Manila Central record 
are due to differences in the averaging times and in 
the exposure, elevation and calibration of the instru- 
ments, as well as to possibly inaccurate estimates of 
the maximum speed in Manila and Pasay City in 1970 
(200 km/hr, see table 2). 

The estimated wind speeds at Baguio based upon the 
1950-73 record are higher than those obtained from 
the 1914-40 data. No explanation is offered for these 
differences; an investigation into their causes seems 
warranted. 

The record at Calapan illustrates the limitations of the 
approach to the definition of design wind speeds 
based on the statistical analysis of the highest annual 
winds. From the data covering the period 1961-72, the 
estimated 50-yr wind based on a Type I distribution is 
88 mph (141 km/hr) [12], versus 1 31 mph (209 km/hr), 



as obtained if the data covering the period 1959-1973 
are used (see table 3). Since wind loads are propor- 
tional to the square of the wind speeds, the ratio be- 
tween the respective estimated winds loads is 
(209/141) 2 = 2.2. 

Although the record at Pasay City is best fitted by a 
type II distribution with opt (y) = 2, it is unlikely, as 
noted previously, that such a distribution correctly 
describes the behavior of the extreme winds. This is 
obvious, particularly in the case of the 1000-yr wind, 
which, on physical grounds, could not possibly attain 
509 mph (820 km/hr) (see table 3). 

The National Structural Code of the Philippines 
specifies, for Zone II and elevations under 9.15 m, a 
design wind of 109 mph (175 km/hr). In the light of 
the data shown in table 2, the value appears to be 
reasonable. It will be noted that tables 2 and 3, and 
figure 13 of Reference 15 indicate that the extreme 
speeds and the frequency of occurrence of tropical 
cyclones, are considerably higher at Laoag than at 
Cebu. This suggests that Zone II could be divided, ac- 
cordingly, into two subzones, with wind load require- 
ments higher in the northern than in the southern 
subzone. 

1.7.3 Zone I. 

As indicated previously, if y(opt) is small , i .e., if the 
differences among maximum wind speeds r>_ \>rded 
in various years are large, the probability dis.i/ - \ 
tions that best fit the data may not describe co; 
the extreme wind speeds for large recurrence im<?r 
vals. The minimum and the maximum winds for the 
period of record are, at Legaspi, 25 mph (40 km/hr) 
and 127 mph (204 km/hr), respectively, and, at 
Tacloban, 26 mph (42 km/hr) and 120 mph (194 
km/hr), respectively. In the writer's opinion, the 
reliability of the N-year wind estimates obtained at 
these stations for N=50, 100 and 1000 is therefore 
doubtful. The same comment applies to the estimates 
for Infanta, where the record length is quite insuffi- 
cient (14 years). The writer therefore believes that the 
results of table 3 should not be used to assess the ade- 
quacy of the design wind speed requirement for Zone 
I specified in Reference 15. Rather, it is reasonable to 
base such an assessment on a comparison between 
wind speeds in Zone I and in areas affected by hur- 
ricanes in the United States. In the light of U.S. ex- 
perience, it is the opinion of the writer that from such 
a comparison it follows that the 124 mph (200 km/hr) 
wind speed requirement for Zone I and elevations 
under 30 ft. (9.15 m) is adequate for structural design 
purposes. 
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1.8 CONCLUSIONS 



ACKNOWLEDGMENTS 



From the analysis of available extreme speed data in 
the Philippines, the following conclusions may be 
drawn: 

1. The design wind speeds specified by the National 
Structural Code of the Philippines for the Zone III 
part of Mindanao appear to be conservative and 
might be somewhat reduced. For this conclusion to 
be validated, it would be necessary to determine, 
from long-term records of tropical cyclone occur- 
rences, that the data analyzed herein are represen- 
tative for southern Mindanao. 

2. Methodological difficulties and uncertainties with 
regard to the reliability of the data preclude, at this 
time, the estimation for Zones II and I of N-year ex- 
treme winds that could be used, with a sufficient 
degree of confidence, as design values within the 

. framework of an explicitly probabilistic code. 

3. According to the data included herein. Zone II can 
be divided into two subzones, with wind load re- 
quirements higher in the northern than in the 
southern subzone. 

4. The data included herein suggest that the wind 
speed requirement specified by the National Struc- 
tural Code of the Philippines for Zone I is adequate 
for purposes of structural design, except as noted 
below. 

5. Higher wind speed values than those specified by 
the National Structural Code of the Philippines 
should be used — except perhaps in the Zone 111 part 
of Mindanao — in open, and in coastal exposure. 

6. Improved design criteria for Zones II and I, includ- 
ing possible redefinitions of these zones, could in 
the future be achieved by applying the 
methodology briefly described at the end of the 
section "Assessment of Procedure Based on the 
Highest Average Monthly Speed." This would re- 
quire, in addition to data on the frequency of oc- 
currence of tropical cyclones at various locations in 
the Philippines, that the following data be availa- 
ble: 

a. Reliable wind speeds, carefully defined with 
respect to terrain roughness, averaging time and 
distance from shore line. 

b. Approximate radial dimensions of tropical 
cyclones. 

c. Approximate dependence of tropical cyclone 
speeds upon minimum central pressure and dis- 
tance from storm center. 
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FIGURE 1. RATIO, r, OF MAXIMUM PROBABLE WIND SPEEDS AVERAGED OVER t SECONDS 
TO THOSE AVERAGED OVER 2 SEC. 
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FIGURE 3b. TYPE I DISTRIBUTION. 
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Facing Page: This wind tunnel at the University of the 
Phillipines is used to study wind effects on scale-model 
buildings. Shown is a model of the CARE, Inc., test house. 
The rows of blocks on the floor of the tunnel generate tur- 
bulence or gustiness similar to that observed in full scale. 
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A GUIDE TO THE 
DETERMINATION OF 
WIND FORCES 
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by simple rules or mathematical formulae. Wide 
variations in building size and shape, type of wind 
exposure, local topography and the random nature of 
wind all tend to complicate the problem. Only by 
direct observation of full-scale situations or by resort- 
ing to properly conducted wind-tunnel tests can the 
characteristics of these flows be established. In spite of 
these complications, guidance can be provided by 
considering some typical flow situations. 




2.2.1 Typical Wind Flow Around Buildings 

A typical flow situation is illustrated in figure 4 
where the wind is blowing face-on to a building with 
a gable roof. The flow slows down or decelerates as it 
approaches the building, creating a positive pressure 
on the windward face. Blockage created by the build- 
ing causes this flow to spill around the corners and 
over the roof: The flow separates (becomes detached 
from the surface of the building) at these points and 
the pressure drops below atmospheric pressure, creat- 
ing a negative pressure or suction on the endwalls and 
on certain portions of the roof. 

A large low-pressure region of retarded flow is cre- 
ated downwind of the building. This region, called 
the wake, creates a suction on the leeward wall and 
leeward side of the roof. The pressures are neither 
uniform nor steady due to the turbulent character of 
the oncoming wind and varying size and shape of the 




FIGURE 4. TYPICAL FLOW PATTERN AND 
SURFACE PRESSURES. 

wake. However, it has been established that the pat- 
terns of wind flow around bluff bodies such as the 
building in figure 4 do not change appreciably with a 
change in wind speed. 

This allows dimensionless pressure coefficients (to be 
discussed later) determined for one wind speed to be 
applied to all wind speeds. In general, the wind 
pressure is a maximum near the center of the wind- 
ward wall and drops off rapidly near the corners. 
Pressures on the side or endwalls are also non- 
uniform, the most intense suctions occurring just 
downstream of the windward corners. 

2.2.2 Effect of Roof Slope 

The pressures acting on a roof are highly dependent 
upon the slope of the roof, generally being positive 
over the windward portion for slopes greater than 30 
degrees. For slopes less than 30 degrees, the wind- 
ward slope can be subjected to severe suctions which 
reach a maximum at a slope of approximately 10 
degrees. Under extreme wind conditions, these suc- 
tions can be of sufficient intensity to overcome the 
dead weight of the building, thus requiring a positive 
tiedown or anchorage system extending from the roof 
to the foundation to prevent loss of the roof system or 
uplift of the entire building. 

Intense suctions are likely to occur along the edges of 
roofs and along ridge lines due to separation or 
detachment of the flow at these points. For certain 
combinations of roof slope and wind direction, a coni- 
cal vortex can be developed along the windward 
edges'of the roof as shown in figure 5. This isa "roll- 
ing up" of the flow into a helical pattern with very 
high speeds and, consequently, very intense suctions. 
If not adequately provided for in the design, these 
vortices along the edges of the roof can cause local 
failures of the roofing, often leading to complete loss 
of the roof. Areas where intense suctions can be ex- 
pected are shown in figure 6. 

2.2.3 Roof Overhangs 

In calculating the total uplift load on a roof, the 
pressure acting on the underside of roof overhangs 
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must also be included. These pressures are usually 
positive and the resultant force acts in the same direc- 
tion as the uplift force due to suction on the top sur- 
face of the roof. Pressures acting on the inside of the 
buiding (to be discussed later) can also contribute to 
the total uplift force and must likewise be accounted 
for. 



VORTICES PRODUCED 
ALONG EDGE OF ROOF. 
WHEN WIND BLOWS 
ON TO A CORNER 



FIGURE 5. VORTICES ALONG EDGE OF 
ROOF. 





AREAS WHERE HIGH 
SUCTIONS MUST BE 
ALLOWED FOR ON THE 
CLADDING 



FIGURE 6. AREAS OF INTENSE SUCTIONS. 



2.3 DESIGN WIND SPEED 

Several factors must be considered in selecting a wind 
speed on which to base the design loads for a building 
or other structure. These include the climatology of 
the geographic area, the general terrain roughness, 
local topographical features, height of the building, 
expected life of the building and acceptable level of 
risk of exceeding the design load. The assessment oi 
climatological wind data and the procedure for ob- 
taining basic wind speeds are discussed in section 1 '). 
The selection of the basic wind speed and the deter- 
mination of modifying factors to obtain the design 
wind speed are discussed in the following sections. 

2.3.1 Mean Recurrence Interval 

The selection of a mean recurrence interval, with 
which there is associated a certain basic wind speed, 
depends upon the intended purpose of a building and 
the consequences of failure. The mean recurrence in- 
tervals in table 4 are recommended for the various 
classes of structures. 

2.3.2 Risk Factor 

There is always a certain risk that wind speeds in ex- 
cess of the basic wind speed will occur during the ex- 



pected life of a building. For example the probability 
that the basic wind speed associated with a 50-year 
mean recurrence interval will be exceeded at least 
once in 50 years is 0.63. The relationship between risk 
of occurrence during the expected building life and 
the mean recurrence interval is given in table 5. It 
should be noted that the risk of exceeding the basic 
wind speed is, in general, not equal to the risk of 
failure. 

2.3.3 Averaging Time and Peak Wind Speed 

It is well known that the longer the time interval over 
which the wind speed is averaged, the lower the indi- 
cated peak wind speed will be. The calculated design 
loads will thus depend upon the averaging time used 
to determine the design wind speeds. In this docu- 
ment, it has been assumed that all speeds used in 
pressure and load calculations are based upon an 
averaging time of 2 seconds. Wind speeds for averag- 
ing times other than 2 seconds can be converted into 
2-second average speeds using the procedure 
described in section 1.0. 

2.4 DESIGN PRESSURES 
2.4.1 Dynamic Pressure 

When a fluid such as air is brought to rest by impact- 
ing on a body, the kinetic energy of the moving air is 
converted to a dynamic pressure q, in accordance 
with the formula 



q=\!2fiU 2 



(1) 



where q = N/m 2 , p is the mass density of the air in 
kg/m 1 and U is the free-stream or undisturbed wind 
speed in m/s. The mass density of air varies with tem- 
perature and barometric pressure, having a value of 
1.225 kg/m 1 at standard atmospheric conditions. In 
the case of tropical storms, the mass density may be 5 
to 10 percent lower. However, this is offset somewhat 
by the effect of heavy rainfall, and the value quoted 
above should be used for all wind pressure calcula- 
tions, i.e., 

<f = 0.613l/ 2 (2) 

2.4.2 Mean and Fluctuating Components of 
Pressure 

As in the case of wind speed, pressures acting on a 
building are not steady, but fluctuate in a random 
manner about some mean value. A typical recording 
of wind speed and pressure at a point on the roof of a 
house is shown in figure 7. 

A close inspection of figure 7 reveals the following 
characteristics: 
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(a) The average or mean pressure is negative (suc- 
tion) 

(b) Pressure fluctuations tend to occur in bursts 

(c) Maximum departures from the mean are in the 
negative (suction) direction 

(d) The peak values far exceed the mean value 




FIGURE 7. TYPICAL RECORD OF WIND 
SPEED AND SURFACE 
PRESSURE. 



To quantify these pressures, it is essential that a suffi- 
ciently long time interval be used to obtain a stable 
mean, p. The fluctuations are described by their stan- 
dard deviation or root-mean-square, p rms . taken 
about the mean. Finally, the peak pressure fluctua- 
tions are described by a peak factor, g, which indi- 
cates the number of standard deviations that the peak 
pressure deviates from the mean. Thus, the peak 
pressure can be expressed as 



J'max =/' + 8 J'rms wr /'mi n = P"« Prms 



(3) 



It should be noted that the peak factor, g, is a random 
variable and has a probability distribution function 
that depends on the geometry of the building and tur- 
bulent structure of the wind. The values of g are 
selected so that the associated probabilities of being 
exceeded are in line with the expected building life. 

2.4.3 Pressure Coefficients 

It is convenient to express pressures acting on the sur- 
faces of a building in terms of the dynamic pressure as 
follows 



pressure, give the wind loads acting in a direction 
normal to those surfaces or elements. The total resul- 
tant forces and moments acting on a building can then 
be determined by considering the appropriate compo- 
nents of these loads acting on each of the surfaces or 
elements. 

As discussed in the previous section, the instan- 
taneous peak pressure can be expressed in terms of the 
mean pressure and a fluctuating component. Since 
pressure fluctuations are limited in spatial extent, it is 
necessary to consider the size of the building surface 
or element when selecting the pressure coefficient. 

Pressures on Extended Areas: For the purpose of 
determining wind loads acting on sizeable surface 
areas such as the walls and roof of a building, the 
pressure coefficients listed in tables 6 and 7 should be 
used. These coefficients have been determined experi- 
mentally from measurements taken on full-scale 
buildings and from wind-tunnel tests and they repre- 
sent an upper limit of conditions likely to occur on the 
indicated building surfaces. 

Pressures on Localized Areas: It is to be expected that 
the smaller the area considered, the larger the effec- 
tive peak pressure will be. In addition, there are cer- 
tain surface areas where intense suctions occur as 
pointed out in sections 2.2.1 and 2.2.2. To provide for 
these cases, pressure coefficients for localized areas 
are included in tables 6 and 7. These coefficients are 
for the purpose of assessing wind loads on local clad- 
ding and roofing elements and should not be used to 
calculate overall loads on buildings. They should be 
used in conjunction with the internal pressure coeffi- 
cients (where appropriate) as described in the follow- 
ing. 

Internal Pressures: As indicated in section 2.2.3 the 
net load or force acting on the roof or walls of a build- 
ing depends not only on the external surface 
pressures, but on the internal pressure as well. The 
magnitude of the internal pressure depends upon the 
building geometry, size and location of openings, and 
wind speed and direction. As with external pressures, 
it is convenient to express internal pressures in terms 
of the dynamic pressure and a pressure coefficient 
C„i- These coefficients can be positive or negative as 
indicated in table 8. The net pressure acting on a 
building element is the algebraic sum of the external 
and internal pressures 



P = Cpfl 



(4) 



where Cp is a pressure coefficient whose value de- 
pends upon the geometry of the building and local 
flow conditions. Pressure coefficients are specified for 
particular surfaces or elements of a building and, 
when multiplied by the surface area and dynamic 



p = q(C p -C pi ) 



(5) 



Thus a positive internal pressure will increase the 
loading on those areas of roofs and walls subjected to 
external suction. 
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2.4.4 Correction Factor for Height of Building 

The pressure coefficients described above are based on 
building heights of 33 ft (10 m) and peak wind speeds 
at 33 ft (10 m) above ground, averaged over 2 seconds. 
Overall loads calculated for buildings appreciably less 
than 33 ft (10 m) in height (measured to eaves or 
parapet) will thus be overestimated if these coeffi- 
cients are used without modification. On the other 
hand, tributary areas such as doors, windows, clad- 
ding and roofing elements will respond to pressure 
fluctuations with duration times considerably less 
than 2 seconds. To account for this, the pressures must 
be multiplied by the correction factors, R, in table 9. 
Thus the expression for the net pressure acting on a 
building surface becomes 

p = <,(C p R-Cp / R / ) (6) 

and the force acting normal to a surface of area A is 

F=(?(CpR-C pi R,-)A (7) 

where R and Rj are correction factors for external and 
internal pressures, respectively. 

2.5 PROCEDURE FOR CALCULATING 
WIND FORCES 

The procedure for calculating wind forces on a build- 
ing is summarized in the following steps. 

1. Select the appropriate mean recurrence interval 



from table 4 

2. Check the associated factor of risk in table 5 and 
select a longer mean recurrence interval if ap- 
propriate. 

3. Determine the basic wind speed for this mean 
recurrence interval and the appropriate terrain 
roughness and type of exposure as outlined in sec- 
tion 1. 

4. Convert the resulting basic wind speed to a 2-se- 
cond mean speed using the procedure described in 
section 1. 

5. Calculate the dynamic pressure q using the ex- 
pression 

17=0.613 U 1 

6. Select the appropriate pressure coefficients from 
tables 6, 7 and 8. 

7 Select the appropriate correction factors from ta- 
ble 9. 

8. Calculate the pressures from the expressions 

P=<?CpR 
or 

p-ftCpR-Cp/Rj) 

9. Multiply these pressures by the respective surface 
areas to obtain the wind forces. 

10. Sum appropriate components of these forces to 
obtain net uplift and drag loads. 
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TABLE 4 MEAN RECURRENCE INTERVAL 



Class of structure 


Mean recurrence interval years 


All structures other than those set out below. 


50 


Structures which have special post-disaster functions, e.g. hospitals, communications build- 
ings, eic. 


100 


Structures presenting a low degree of hazard to life and other property in the case of 
failure. 


20 



TABLE 5. RELATIONSHIPS BETWEEN RISK OF OCCURRENCE, MEAN RECURRENCE INTER- 
VAL AND EXPECTED LIFE OF BUILDING 



Desired 




Risk of exceeding in N years the wind speed corresponding 




Lifetime 




to the indicated mean recurrence interval 




N 

Years 


0.632 


0.50 


0.40 


0.30 


0.20 


0.10 




Mean Recurrence Interval in Years 


10 


10 


15 


20 


29 


45 


95 


20 


20 


29 


39 


56 


90 


190 


50 


50 


72 


98 


140 


224 


475 


100 


100 


144 


196 


280 


448 


949 



Note: From this table it will be seen that there isa 10'^ risk that the wind speed corresponding to a mean recurrence interval of 475 
years will be exceeded in a lifetime of 50 years. 
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TABLE 6. PRESSURE COEFFICIENTS FOR WALLS OF RECTANGULAR BUILDINGS 



RuiMino 

UU11U 1.1 111 

Height/Width 
Ratio 


DuiiaiiiK 
Length/width 
Ratio 


Winn Anola 

tt inu r\ngit? 

a 

(Degrees) 


Cp for Face 


Local 

C_ 
P 


A 


B 


C 


D 









0.8 


-U.5 


-0.6 


-0.6 






1 <//w<I.5 












-1.2 






911 


-0.6 


-0.6 


0.8 


-0.5 




h/w<0.5 























0.8 


-0.5 


-0.8 


-0.8 






l.5<//w<4 












-1.3 






90 


-0.4 


-0.4 


0.8 


-0.4 








I) 


O.S 


-O.S 


-0.7 


-0.7 






1 <//w<1.5 












-1.4 






un 


-0.7 


-0.7 


0.8 


-0.5 




0.5<h/w<1.5 























0.8 


-0.6 


-0 9 


-0.9 






1.5<//w<4 












-1.4 






90 


-0.4 


-0.4 


0.8 


-0.3 











0.8 


-0.5 


-0 8 


-0.8 






l<//w<1.5 












-1.5 






90 


-0.8 


-0.8 


0.8 


-0.5 




1.5<h/w<4 























0.8 


-0.6 


-0.9 


-0.9 






1.5<//w<4 












-1.5 






90 


-0.4 


-0.4 


0.8 


-0.3 





Notes: (l)h is the height to eaves or parapet, / is the greater plan dimension of the building and w is the lesser plan dimension. 
(2) Local C_ values (last column) should be used in conjunction with correction factors for overall areas in Table 9. 
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TABLE 7. PRESSURE COEFFICIENTS FOR ROOFS OF RECTANGULAR BUILDINGS 



Building 


Wind Angle 


Area 




Roof Slope 




Height/Width 
Ratio 


a 


Designation 




Degrees 




(Degrees) 





10 


20 


25 






EF 


-1.0 


-1.0 


-0.4 


-0.3 






GH 


-0.6 


-0.6 


-0.8 


-0.6 







J 


-1.6 


-1.9 


-1.9 


-!.6 






K 


-1.4 


-1.4 


-2.0 


-1.6 


h/w<0.5 




























EG 


-1.2 


-1.1 


-1.0 


-0.8 






FH 


-0.6 


-0.6 


-0.6 


-0.5 




90 


J 


-2.0 


-1.8 


-1.8 


-1.6 






K 


-1.4 


-1.4 


-2.0 


-1.6 






EF 


-1.0 


-1.0 


-0.6 


-0.4 






GH 


-0.6 


-0.6 


-0.8 


-0.8 







J 


-1.8 


-1.8 


-1.8 


-1.8 


0.5<h/w < 4 




K 


-1.4 


-1.4 


-2.0 


-1.6 












-0.8 




EG 


-1.2 


-1.1 


-1.0 






FH 


-0.6 


-0.6 


-0.5 


-0.4 




90 


) 


-1.8 


-1.8 


-i.6 


-1.6 






K 


-1.4 


-1.6 


-1.6 


-1.6 



Notes: (1) The pressure coefficient on the underside of roof overhangs should be taken as that on the adjoining wall surface. 



(2) Local C v values (/ and K) should be used in conjunction with correction factors for overall areas in table 9. 




Extended areas Local areas 



ELEVATION PLAN 
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TABLE 8. INTERNAL PRESSURE COEFFICIENTS FOR RECTANGULAR BUILDINGS 



Condition 


Internal pressure coefficient Cpj 


1. Two opposite walls equally permeable, other walls imper- 
meable: 

(a) Wind normal to permeable wall 

(b) Wind normal to impermeable wall 


+0.3 
-0.3 


2. Four walls equally permeable 


-0.3 or +0.2 whichever is the more severe for 
combined loadings 


}. Dominant opening on one wall, other walls of equal per- 
meability: 

(a) Dominant opening on windward wall, having a ratio 
of permeability of windward wall to total permeability of 
other walls and roofs subject to external suction, equal 
to— 

2 

3 

f> or more 

(b) Dominant opening >n leeward wall 

(c) Dominant opening on side wall 

(d) Dominant opening in a roof segment 


+0.5 
+0.6 
+0.8 

value of Cj, for leeward external wall surface 

value of C„ for side external wall surface 
value of Cp for external surface of roof segment 



Notes: ( I ) Internal pressures developed within an enclosed structure may be positive or negative depending on the position and 
size of the openings. 



(2) In the context of table 8 the permeability of a surface is measured by the total area of openings in the surface under con- 
sideration. 

(3) The value of C„ ( can be limited or controlled to advantage by deliberate distribution of permeability in the wall or roof, 
or by the deliberate provision of a venting device which can serve as a dominant opening at a position having a suitable 
external pressure coefficient. An example of such is a ridge ventilator on a low-pitch roof, and this, under all directions 
of wind, can reduce the uplift force on the roof. 



TABLE 9. CORRECTION FACTOR (R) FOR HEIGHT OF BUILDING 



Terrain 


Structural System 


Area 


h < 5 


5 < h < 10 


Smooth 

Zo < 0.12 m 


Walls 


Overall 
Elements 


0.85 
1.00 


1.00 
1.20 


Roofs 


Overall 
Elements 


0.85 
1.05 


1.00 
1.25 




Internal Pressure 




0.85 


1.00 


Rough 


Walls 


Overall 
Elements 


0.75 
0.90 


1.00 
1.20 


Zo > 0.12 m 


Roofs 


Overall 
Elements 


0.75 
0.95 


1.00 
1.25 




Internal Pressure 




0.75 


1.00 



Notes: (l)The term "Overall" refers to the entire area of a given wall or roof slope. 



(2) The term "Elements" refers to roof and cladding elements, doors, windows, etc. 

(3) The terrain roughness parameter Z must be estimated subjectively. The following values are suggested for various 
types of exposure. 

TYPE OF EXPOSURE Z„ (meters) 

Coastal 0.005-0.01 

Open country 0.02-0.12 

Outskirts of towns, suburbs 0.13-*3.3u 

Centers of towns 0.40 



APPENDIX A 

ILLUSTRATIVE EXAMPLE 

A housing development is to be located in flat, open 
country on the outskirts of Zamboanga, Philippines, 
and will ultimately consist of several hundred single- 
family dwellings of quite similar geometry. The 
period of construction is anticipated to be from 10 to 
15 years. The basic plan dimensions are 6.2 x 7.5 m 
and the height to the eaves is 2.7 m. The gable roof has 
an overhang of 0.7 m on all sides and a slope of 1 
degrees. Openings for doors and windows are evenly 
distributed on the exterior walls. 

Because the development is to be built over a period of 
several years, it would not be appropriate to assume a 
built-up area in selecting the basic wind speed and 
flat, open country will be assumed here. 

From table 4, a mean recurrence interval of 50 years 
is selected and it is considered that the associated risk 
of exceeding the basic wind speed (0.632) in table 5 is 
acceptable. 

From section 1, the 1-minute average wind speed 
(N=50) for Zamboanga is 88 km/hr (Type I distribu- 
tion). Since this is based on data obtained in open 
country at 10 m above ground, this speed can be con- 
verted directly to the design speed. Also, from sec- 
tion 1 the ratio of the 1-minute speed to the 2-second 
peak speed is 0.82. Thus the design speed is 

U = 88/0.82 = 107.3 km/hr = 29.8 m/s 

The dynamic presure is calculated from equation 2 of 
section 2.4.1 

q = 0.613U 2 = (0.613) (29.8)* = 544 N/m 2 

Wind pressures are next calculated using equations 
4-6 and the coefficients presented in table 6-9. Note 
that 

h/w = 2.7/6.2 = 0.44 

and 

Z/w = 7.5/6.2 = 1.21 

WALLS 

Inspection of tables 6 and 8 reveals that the worst 
cases are walls A and C with the wind blowing nor- 
mal to the ridge. For wall A, C„ = 0.8 and for wall C, 
Cp = -0.6. The local Cp is -1.2. The internal pressure 
coefficients can range from 0.2 to -0.3. Table 9 indi- 
cates that the reduction factor is 0.85 for walls and in- 
ternal pressures and 1.00 for cladding elements, doors, 
windows, etc. 



For wall A, 

p=(544)l0.8-(-0.3)](0.85) 
= 509 N/m 2 

For wall C, 

p=(544)i-0.6-(0.2)](0.85) 
= -370 N/m 2 

For cladding elements, the worst cases are 

p= (544) 10.8 -(-0.3) (0.85)] 
= 574 N/m 2 



and 

p = (544) 1-0.6 -(0.2) (0.85)1 
= -419 N/m 2 

For local pressures acting on strips of width 0.2 w = 
1 .2 m at each corner, 

p = (544) (-1.2 - 0.2) (0.85) 
= -647 N/m 2 

ROOF 

inspection of table 7 reveals that the greatest uplift 
pressures on extended areas occur when the wind is 
blowing along the ridge. 

For sections E and C, 

p = (544) [-1.1 -(0.2)] (0.85) 
= -601 N/m 2 

For sections F and H, 

p = (544) 1-0.6 -(0.2)] (0.85) 
= -370 N/m 2 

Pressures acting on roofing elements in sections E and 
G a*e obtained as follows: 

p = (544) [(- 1 . 1) ( 1 .05) - (0.2) (0.85)] 
= -727 N/m 2 

and for sections F and H, 

p = (544) K-0.6Ml.05) - (0.2X0.85)] 
= -438 N/m 2 

Localized pressures act on strips of width 0.15 w = 
0.93 m as shown in table 7. The worst case occurs for 
area J with the wind blowing normal to the ridge. 
Note that the uplift pressure under the eaves must 
also be included. 
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p = (544) 1-1.9 -(0.8)1 (0.85) 
= -1.2k N/m 1 



COMMENT 



For area K in section F, this negative pressure or suc- 
tion is slightly less 

p = (544) 1-1.4 -(0.8)] (0.85) 
= -1.0* N/m* 

Along the ridge (area K), the localized pressure is 

p = (544) [-1.4 (0.2)] (0.85) 
= -740 N/m 3 

TOTAL UPLIFT FORCE 

The total uplift force on the building is calculated for 
the wind blowing normal to the ridge as follows: 

Area of one roof slope = [7.5 +(2) (0.7)] [6.2/ (2Cos 10°) 

+ 0.7] 
= (8.9)0.85) 
= 34.2 m 2 

Note that areas E, F, G and H include areas J and K 
when calculating overall loads. 

Uplift = (544) (1.0 + 0.6) (34.2) (Cos 10 Q ) (0.85) 
+ (544) (6.2) (7.5) (0.2) (0.85) 

= 29.2kN 
TOTAL DRAG FORCE 



The loads calculated above are the loads that can 
reasonably be expected to occur under the conditions 
stated in the example. They should be considered as 
the minimum suitable loads for use with stresses and 
load factors appropriate for the type of structural 
material used. 

For geographical areas exhibiting large variations in 
annual extreme wind speeds, the basic wind speed 
should be selected with caution. The application of 
probabilistic models of extreme wind speeds and 
some of their limitations are discussed in section 1.0. 



The total drag force (neglecting the roof) is calculated 
as the sum of the loads on the windward and leeward 
walls. 

Drag = (544) (2.7) (7.5) 10.8 - (-0.5)] (0.85) 
= 12.2kN 
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ANNOUNCEMENT OF NEW PUBLICATIONS IN 
BUILDING SCIENCE SERIES 

Superintendent of Documents, 
Government Printing Office, 
Washington, D.C., 20402 

Dear Sir: 

Please add my name to the announcement list of new publications to be 
issued in the series: National Bureau of Standards Building Science Series. 

Name - - 

Company 

Address - 

City State Zip Code 

(Notification key N-339) 



NBS TECHNICAL PUBLICATIONS 



PERIODICALS 

JOURNAL OF RESEARCH reports National Bureau 
of Standards research and development in physics, 
mathematics, and chemistry- It is published in two 
sections, available separately: 

• Physics and Chemistry (Section A) 

Papers of interest primarily to scientists working in 
these fields. This section covers a broad range of physi- 
cal and chemical research, with major emphasis on 
standards of physical measurement, fundamental con- 
stants, and properties of matter. Issued six times a year. 
Annual subscription: Domestic, $17.00; Foreign, $21.25. 

• Mathematical Sciences (Section B) 

Studies and compilations designed mainly for the math- 
ematician and theoretical physicist. Topics in mathemat- 
ical statistics, theory of experiment design, numerical 
analysis, theoretical physics and chemistry, logical de- 
sign and programming of computers and computer sys- 
tems. Short numerical tables. Issued quarterly. Annual 
subscription: Domestic, $9.00; Foreign, $11.25. 

DIMENSIONS/NBS (formerly Technical News Bulle- 
tin) — This monthly magazine is published to inform 
scientists, engineers, businessmen, industry, teachers, 
students, and consumers of the latest advances in 
science and technology, with primary emphasis on the 
work at NBS. The magazine highlights and reviews 
such issues as energy research, fire protection, building 
technology, metric conversion, pollution abatement, 
health and safety, and consumer product performance. 
In addition, it reports the results of Bureau programs 
in measurement standards and techniques, properties of 
matter and materials, engineering standards and serv- 
ices, instrumentation, and automatic data processing. 

Annual subscription: Domestic, $12.50; Foreign,$15.65. 

N0NPERI0DICALS 

Monographs — Major contributions to the technical liter- 
ature on various subjects related to the Bureau's scien- 
tific and technical activities. 

Handbooks — Recommended codes of engineering and 
industrial practice (including safety codes) developed 
in cooperation with interested industries, professional 
organizations, and regulatory bodies. 
Special Publications — Include proceedings of conferences 
sponsored by NBS, NBS annual reports, and other 
special publications appropriate to this grouping such 
as wall charts, pocket cards, and bibliographies. 
Applied Mathematics Series — Mathematical tables, man- 
uals, and studies of special interest to physicists, engi- 
neers, chemists, biologists, mathematicians, com- 
puter programmers, and others engaged in scientific 
and technical work. 

National Standard Reference Data Series — Provides 
quantitative data on the physical and chemical proper- 
ties of materials, compiled from the world's literature 
and critically evaluated. Developed under a world-wide 
program coordinated by NBS. Program under authority 
of National Standard Data Act (Public Law 90-396). 



NOTE: At present the principal publication outlet for 
these data is the Journal of Physical and Chemical 
Reference Data (JPCRD) published quarterly for NBS 
by the American Chemical Society (ACS) and the Amer- 
ican Institute of Physics (AIP). Subscriptions, reprints, 
and supplements available from ACS, 1155 Sixteenth 
St. N.W., Wash. D. C. 20056. 

Building Science Series — Disseminates technical infor- 
mation developed at the Bureau on building materials, 
components, systems, and whole structures. The series 
presents research results, test methods, and perform- 
ance criteria related to the structural and environmental 
functions and the durability and safety characteristics 
of building elements and systems. 
Technical Notes — Studies or reports which are complete 
in themselves but restrictive in their treatment of a 
subject. Analogous to monographs but not so compre- 
hensive in scope or definitive in treatment of the sub- 
ject area. Often serve as a vehicle for final reports of 
work performed at NBS under the sponsorship of other 
government agencies. 

Voluntary Product Standards — Developed under proce- 
dures published by the Department of Commerce in Part 
10, Title 15, of the Code of Federal Regulations. The 
purpose of the standards is to establish nationally rec- 
ognized requirements for products, and to provide all 
concerned interests with a basis for common under- 
standing of the characteristics of the products. NBS 
administers this program as a supplement to the activi- 
ties of the private sector standardizing organizations. 
Consumer Information Series — Practical information, 
based on NBS research and experience, covering areas 
of interest to the consumer. Easily understandable lang- 
uage and illustrations provide useful background knowl- 
edge for shopping in today's technological marketplace. 

Order above NBS publications from: Superintendent 
of Documents, Government Printing Office, Washington, 
D.C. 20402. 

Order following NBS publications— NBSIR's and FIPS 
from the National Technical Information Services, 
Springfield, Va. 22161. 

Federal Information Processing Standards Publications 
(FIPS PUBS)— Publications in this series collectively 
constitute the Federal Information Processing Stand- 
ards Register. Register serves as the official source of 
information in the Federal Government regarding stand- 
ards issued by NBS pursuant to the Federal Property 
and Administrative Services Act of 1949 as amended, 
Public Law 89-306 (79 Stat. 1127), and as implemented 
by Executive Order 11717 (38 FR 12315, dated May 11, 
1973) and Part 6 of Title 15 CFR (Code of Federal 
Regulations). 

NBS Interagency Reports (NBSIR) — A special series of 
interim or final reports on work performed by NBS for 
outside sponsors (both government and non-govern- 
ment). In general, initial distribution is handled by the 
sponsor; public distribution is by the National Techni- 
cal Information Services (Springfield, Va. 22161) in 
paper copy or microfiche form. 



BIBLIOGRAPHIC SUBSCRIPTION SERVICES 



The following current-awareness and literature-survey 
bibliographies are issued periodically by the Bureau: 
Cryogenic Data Center Current Awareness Service. A 
literature survey issued biweekly. Annual subscrip- 
tion: Domestic, $25.00 ; Foreign, $30.00). 
Liquified Natural Gas. A literature survey issued quar- 
terly. Annual subscription: $20.00. 



Superconducting Devices and Materials. A literature 
survey issued quarterly. Annual subscription: $30.00'. 
Send subscription orders and remittances for the pre- 
ceding bibliographic services to National Bureau of 
Standards, Cryogenic Data Center (275.02) Boulder, 
Colorado 80302. 



U.S. OCPMtTMENT OP CO MM E nC C 
National Bmtmh of I 

Washington. O.C. 80234 



OFFICIAL BUSINESS 

Pinalty lor Private Uu. S300 



POSTASE AND FEES PAIO 
U.S. DEPARTMENT OF COMMERCE 
COM-21B 



SPECIAL FOURTH-CLASS RATE 
BOOK 




